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COUPLING PATTERN OF VARIABLES 



JEAN B. LASSERRE 

Abstract. We specialize Schmiidgen's Positivstellensatz and its Putinar and 
Jacobi and Prestel refinement, to the case of a polynomial / € M.[X, Y] + 
lR[y,Z], positive on a compact basic semi-algebraic set K described by poly- 
nomials in R[X, Y] and R[Y, Z] only, or in R[X] and R[Y, Z] only (i.e. K 
is cartesian product). In particular, we show that the preordering P{g,h) 
(resp. quadratic module Q{g, h)) generated by the polynomials {gj} C R[X, Y] 
and {hj^} C lR[y, Z] that describe K, is replaced with P(g) + P{h) (resp. 
Q{g)+Q{h)), so that the absence of coupling between X and Z is also preserved 
in the representation. A similar result applies with Krivine's Positivstellensatz 
involving the cone generated by {gj,hk}. 



1. Introduction 

Schmudgen's Positivstellensatz 8 and its Putinar [7j and Jacobi and Prestel 

refinements, are very useful tools with (relatively recent) particular impact in 
polynomial optimization; see e.g. Lasserre 0, Schweighofer j^. 

An interesting issue, and not only from a computational viewpoint, is to derive 
a Positivstellensatz that preserves a possible coupling pattern of variables present 
in the original polynomial / (positive on a compact semi-algebraic set). By this 
we mean that if there is no coupling of variables Xi and Xj in / as well as in the 
polynomials that describe the basic semi-algebraic set, we would like to obtain a 
representation in which the same property holds. 

A first positive result in this vein is derived in Lasserre ^ under some condition 
of the coupling pattern, known as the running intersection property in graph theory. 
Such specialized representations are particularly important from a computational 
viewpoint, as evidenced by the impressive computational experiments presented in 
Waki et al. when used in polynomial optimization problems with structured 
sparsity. 

Here we present two specialized Positivstellensatz when / S M.[X, Y] + M[Y, Z], 
that is, when there is no coupling between variables X and Z in /. While in the 
first one, the compact basic semi-algebraic set K is also described by polynomials in 
R[X, Y] and M.[Y, Z] only, in the second one K is described by polynomials in M.[X] 
and R[Y, Z] only, i.e., K — K.^ x for some K^;, K^^. Our result does not require 
any assumption at all (except of course compactness of K). Although part of our 
result in the first case, namely Theorem l2.1f b^ below, could be derived from with 
appropriate modifications, the general form of our specialized Positivstellensatz in 
Theorem l2.1f a'l-(b')-(c') is not apparent from ^6^, and we think it is important enough 
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and of self-interest to deserve a special treatment and presentation to an audience 
not necessarily aware of the more computational oriented result _6 . On the other 
hand, the more involved case where K is a cartesian product cannot be deduced 
from 121 . 

Contribution. Let K^y C M"+™, Ky^ C M™+p, and K C M'^+^+p be basic 
compact semi-algebraic sets defined by 

(1.1) K,y 
(1-2) Ky, 

(1.3) K 

for some polynomials {gj} C R[X, F], {hk} C M[y, Z], and some finite index sets 

'^xy ; Iy2 ^ 

Let P{g) C r] and P{h) C R[Y, Z] be the preordering generated by 

{9j}jei^y and {hk}k(^iy^, respectively. 

Let / € M.[X, Y] + M.[Y, Z], i.e., there is no coupling of variables X and Z in /. 

• We first obtain the following specialized Positivstellensatz. 

(1.4) [/e R[X,r]+M[r,Z] and/ > OonK] => f e P{g) + Pih), 

to compare with Schmiidgen's Positivstellensatz which states that / e P{g, h). 

The Positivstellensatz H1.4II is a specialization of Schmiidgen's Positivstellensatz 
where the preordering P{g,h) is replaced with P{g) + P{h). And so, only poly- 
nomials in M[X, F] and M[y, Z] are involved in the representation H1.4I) . In other 
words, the absence of coupling between the variables X and Z is preserved in the 
Positivstellensatz. 

If in addition to be compact, 'K.xy and/or Kj,^ satisfy Putinar's condition in 
[7], then in 11.4|l one may replace P{g) and/or P{h) with the quadratic modules 
Q{g) C R[X,y] and Q{h) C R[F, Z] generated by {5j}jei,.„ and {hk}k£iy,, respec- 
tively. 

Finally, assume that < < 1 on Kj;j, for all j € Ixy , and < /i^ < 1 on Kj,^ for 
all k € iyz)- If the famihes {0, 1, {gj}} and {0, 1, {hk}} generate the algebra R[X, Y] 
and M.[Y, Z] respectively, then one also obtains the alternative representation 

(1.5) [/ G R[X,Y]+R[Y,Z] and / > on K] ^ / e C(g, 1 - .g) + C(/i, 1 - /i), 

where C'{g, 1 — g) (resp. C(/i, 1 — h)) is the cone generated by the polynomials 
1 - 9j} (resp. {hk, 1 - hk}). 
The Positivstellensatz H1.5|l is a specialization of Krivine ^ and Vasilescu [1(1) 
Positivstellensatz where the cone C{g, h,l — g,l — h) is replaced with C{g, 1 — g) + 
C{h, 1 — h). And so, only polynomials in R[X, Y] and M.[Y, Z] are involved in the 
representation l|1.5|l . In other words, the absence of coupling between the variables 
X and Z is also preserved in the Positivstellensatz H1.5|l . 

• When K^j, is now replaced with K.^ C M", that is, K = K.^ x ^yz, then with 
/ e M.[X, Y] + R[Y, Z], we now obtain 

(1.6) [/ > OonK, xK,,] ^ / e S^^ + P(.g) + P(/i), 

where now P{g) C M.[X] and S^y C M.[X, Y] denotes the set of sums of squares. 
This case is more involved and cannot be derived from 0. 



= {{x,y) eR"+™| gjix,y)>0, j G I,,} 
= {{y,z) eM™+f I hk{y,z)>0, kelyz} 
= { {x, y, z) e M"+"+P I {x, v) E Kxy-, {y, z) e K,, } 
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The paper is organized as follows. Our two results are stated in the next section, 
and for clarity of exposition, their proofs are postponed to section iJ21 

2. Main result 

Let R[X, Y, Z] denote the ring of real polynomial in the variables {Xi, . . . , X„), 
(Yi, . . . , Yjn) and {Zi, . . . , Zp). Let ||a;|| denote the euclidean norm of x G M". 

Let c R[X,Y] (resp. S^^ C R[Y, Z]) be the space of elements of R[X,Y] 
(resp. M.[Y, Z]) that are sum of squares (in short s.o.s.). 

Given a family {ffjljei^^ £ F] (resp. {hk}keiy^ C R[F, Z]) for some finite 
index set Ixy C N (resp. lyz C N), denote by P{g) (resp. P{h)) the preordering 
generated by {gj} (resp. {hk})- That is, a G P{g) if 

(2.1) (7= <Tjgj{= J2 '^jJlaj) with ajeSly VJCl.j^ 

(with the convention gj := Y[jej9j = 1 if >/ = 0), and same thing for P{h) with 
obvious adjustments. Similarly, denote by Q{g) C M.[X, Y] the quadratic module 
generated by {gjjjei^y C R[X, F]. That is, a G Q{g) if 

(2.2) a = ao+Y a-i gj with aj G S^y V j G l^y U {0}, 

and same thing for P{h) with obvious adjustments. 

Finally, denote by C{g,l — g) (resp. C{h,l — h)) the cone generated by the 
family {g^, 1 - g^} (resp.{ft.fc, 1 - hk}), i.e. a G C{g, 1 - g) if 

(2.3) a= ^ c^pg'^il-gf^ ^ H 3? " ff^)* 
for some nonnegative scalar coefficients {cafs}- 

Theorem 2.1. iei K^y C R", K;,^ C M", and Kc M"+™+p 6e i/ie &asic compact 
semi- algebraic sets defined in Hl.lp - fl.ci}} . and assume that K has nonempty inte- 
rior. Let / G K[X, F] + M[y, Z] . 

(a) // / is positive on K then f G -P(.g) + P{h). 

(b) //iV - ||(X,F)f G Qig) and/or N - ||(F,^)|P G Q(/i) /or some sca/ar N, 
and if f is positive on K, then in (a) one may replace P{g) with Q{g) and/or P{h) 
with Q{h). 

(c) Assume that the gj 's and hk 's are normalized, i.e., < gj < 1 on K^^j, for 
all j G Ixy, and < hk < 1 on 'K.yz for all k G lyz- In addition, assume that the 
family {0,1, gj,j G T-xy} (resp. {0,1, hk,k G lyz}) generates the algebra R[X,Y] 
(resp. R[Y, Z]). If f is positive on K then f G C{g, 1 ~ g) + C{h, 1-h). 

The proof is postponed to ii3.2l One may see that in the various representations 
of / in TheoremO;a)-(b)-(c), all polynomials in P{g), P{h), Q{g), Q{h), C{g,l-g) 
and C{h,l — h) are either in F] or M[F, Z], and so, there is no coupling of 
variables X and Z, as in /. In other words, the coupling pattern of variables is 
preserved in each representation. 
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As already mentioned, when N -\\{X,Y)\\^ e Q{g) and N - \\{Y, Z)\\^ G Qih) for 
some scalar N, then Theorem l2.ir b) can be derived from a result stated in [Hj, with 
appropriate modifications. This is because the condition (f .3) in (known as the 
running intersection property in graph theory) is satisfied, by taking Ii := {X, Y), 
h^{Y,Z). 

We next consider the more involved case where K is the cartesian product x 
Ky^, with 

(2.4) K, = { a; e M" : gj{x) > 0, j G I, }, 

for some polynomials {gj} C ]R.[X]. Now, both sets P{g) and Q{g) are subsets of 
M.[X]. 

Theorem 2.2. LefKy^ C M^+p andJix C M" be compact basic semi- algebraic sets 
as defined in fj.gj) and respectively. Assume that K := x ^yz C MJ^+^^+p 

has nonempty interior, and let f S M.[X, Y] + R[Y, Z\. 

(a) // / is positive on K then f E S^y + P{9) + P{h)- 

(b) IfN~\\X\\^ e Qig) and/or N - \\{Y, Z)f e Q{h) for some scalar N, and if 
f is positive on K, then in (a) one may replace P{g) with Q{g) and/or P{h) with 
Qih). 

For a proof see ii3.3l So if we compare Theorem 12 . iT a V fb) and Theorem 12. 2r aV 
(b) we see that in the latter the preordering P{g) and the quadratic module Q{g) 
are now in M.[X] instead of M.[X, Y] in the former. On the other hand we need some 
term of S^y. Finally, Theorem 12. 2r b') cannot be deduced from 

3. Proofs 

We first need to introduce some additonal notation and definitions. 

3.1. Notation and definitions. Let u = (ua/j-y) be a sequence indexed in the 
canonical basis {X'^Y'^ Z'') of F, Z], and define the linear functional : 
R[X, Y,Z]^R to be: 

/(= E fo.p-^X'^Yf' Z^) ^ LM)= E U0fUo.f3j. 

Moment matrix. With a sequence u — (uap-y) indexed in the canonical basis of 
R[X, Y, Z] is associated the moment matrix Mr{u) of order r, and defined by 

Mr(u)((a,(3,j),{a',(3\j')) = Lu{X''+''' Y''+''' Z'>+''' ) - u^+o.',p+p'..^+Y, 

for all a, a' G N",/3,/3' G N",7,7' G W, with \a + f3 + j\ < r, and \a' + f3' + "f'\ < r. 
An infinite sequence u = (uafj-y) has a representing measure if 

Ua;37 = J X'^Y'^Z'' d^i{X,Y,Z), a G N", /3 G N", 7 G NP. 

Let Rri-'^, Y, Z] C R[X, Y, Z] be the space of polynomials of degree at most r, and 
denote by s(r) its dimension. If f = (/„^^) G M"^'') denotes the 
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vector of coeificients of an arbitrary polynomial / G M,- [X, Y, Z] then 

(3.1) (f,M,(M)f) = L„(/2) = I fdfi>0, 

and so, as / S R,. [X, Y, Z] was arbitrary, Mr{u) >z 0. 

We next denote by Mr{u,xy) the moment submatrix obtained from Mr{u) by 
retaining only those rows and columns (a,/?, 7) with 7 = 0. Similarly, denote 
by Mr{u,yz) the moment submatrix obtained from M^iu) by retaining only those 
columns and rows (a, /3, 7) with a = 0. Introducing the subsequences u^^ := (uapo) 
and it^^ :— (uoq/j), notice that Mr{u,xy) is just the moment matrix Mr{u^^) of 
the sequence u^^ indexed in the canonical basis (X^Y^) of F]. Similarly, 
Mr{u, yz) is the moment matrix Mr(u*'^) of the sequence m''^ indexed in the canon- 
ical basis {Y'^Z^) of M.[Y, Z\. 

Localizing matrix. Let g G M[X, Y, Z] and let u ~ (ua/37) be a sequence indexed 
in the canonical basis of M.[X, Y, Z\. The localizing matrix Mr{gu) associated with 
g (= J2abc 9abcX°'Y'' Z'^) and u, is obtained from Mr{u) by 



M,(g^.)((«,/3,7),(a',/5',7')) = L^{g{X,Y, Z) X^+^^' Y^+P' Z^^''' ) 

for all a, a' G N", /3,/3' G N" and 7, 7' G with |a+/?+7| < r and \a'+P'+i\ < r. 
Again, if u has a representing measure /i then 

(3.2) (f,M,(5^/)f) = Lu{gf) - ^ s/'^A^, V/ G F, Z]. 

If fi has its support in the level set {{x,y,z) G W^+"^+p | g[x,y,z) > 0} then 
Mr[gu) > 0. 

As for the moment matrix, one may also define the localizing matrices Mr{g u, xy) 
and Mr{hu,yz) associated with u and g G F], h G K[F, Z], respectively. 

They are obtained from Mr{gu) (resp. Mj-ihu)) by retaining only those rows and 
columns (a,/?, 7) with 7 = (resp. a = 0). They can also be considered as the 
localizing matrix Mr{gu^^) (resp. Mr{hu^'^)) associated with g (resp. h) and the 
subsequence (resp. u^^). 

3.2. Proof of Theorem 12.11 (a) For every J C I^-j,, and depending on parity, let 
2r{gj) or 2r{gj) — 1 be the degree of gj (and similarly for 2r{hK) or 2r{hK) ~ 1). 
Let 2ro := max[deg/, 2 max[r'((7,/), r(/ii^)]], and for r > ro, consider the following 

optimization problem: 



{niin Luif) 
s.t. Mr^rigj){g.jU,Xy) "^0, V J C I^y 
Mr-r(hk) {hK u,yz)yO, 'i K <Zly^_ 
"0 = 1 

where g,j := HjeJ 9j for every J C I^y and Hk := Ilfce/f ^oi" "^very if C 1^^. Re- 
call that 50 = 1 and so, Mr{g%u^xy) is just the moment matrix Mr{u,xy) defined 
in tl3.1l the same holds true for /ig. Denote by inf the optimal value of Qr. 
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Qr is a convex optimization problem called a semidefinite programming problem. 
Up to arbitrary fixed precision, it can be solved in time polynomial in the input 
size of the problem data, and efficient specialized software packages are available. 
For more details the interested reader is referred to e.g. Vandenberghe and Boyd 
|12|. The dual problem Q* of Qr is also a semidefinite program, which reads 



(3.4) q: 



max A 

s.t. 



/ - A = ^ O-J .9.7 + ^ i'KhK 

{aj} C S^y, {^^} c S^^, 

degajgj, degipK^K <2r V J I^y, K C ly^ 
Equivalently, Q* reads 

(3.5) Q; : max { A : / - A e Pr{g) + Pr{h) }, 

where Pr{g) C P{g) denote the set of elements a G P{g) that can be written as 
in H2.1|l and where in addition, deg [aj gj) < 2r for all J C I^.^,, and similarly for 

Pr{h). 

The first important step (i) is to prove that inf Q^. "f /* as r ^ oo, with /* :— 
min { f{x, y, z) \ (x, y, z) G K }, or equivalently, 

(3.6) /* = min { [ f d^: \ \ K) = 0; Vj(K) = 1}, 

V' J 

where the infimum is taken over all Borel probability measures ip on R"+™+p. 

The second important step (ii) is to prove absence of a duality gap between Q^. 
and its dual Q*. The final step (iii) easily follows from (a) and (b). 

Step (i). We first prove that inf < /* for all r > ro. Let {xo, yo, zo) £ K, and 
let n := S(^xo,yo,zo) be the Dirac probability measure at (xo, yo, zo). Let u := {uaf3-y} 
be its (well defined) sequence of moments. Then, obviously, Mr{u,xy) ^ 0, and 
Mr{u, yz) ^ 0. Next, as /i is supported on K, we obviously have Mr{gj u, xy) >r 0, 
for every J C and r £ N. Similarly, Mr{hK u,yz) >: 0, for every K C Ij,^ and 
r G N. Therefore, u is feasible for with value Ly{f) = J fdfi — f{xo,yo,zo), 
which proves that 

(3.7) inf Qr < f* := min { /(x, y, z) \ (x, y,z) eK }, Vr > tq. 

We next prove that inf > ~oo. Let r > rg be fixed. As K^^^ is compact, there 
is some Nr such that the polynomials Nr ± X'^Y'^ , |a + /3| < 2r, are all positive 
on Kxy By Schmiidgen's Positivstellensatz they all belong to the preordering 
P{g). There is even some l{r) such that they all belong to Pi(^j.){g)- Similarly, as 
is compact, there is some t{r) such that the polynomials N'^ — belong 
to Pt(r){h), for some N'^. and all [a + /3| < 2r. 

So, let be an arbitrary feasible solution of Q^, with s > max[/(r), t(r)]. One 
has 

Lu^ (Nr ± X"y^) ^ Nru'o± u%pQ > 0, V |a + /3| < 2r, 
because Nr ± X°'Y^ = X],7ci^„ ^j9j with degcrj gj < 2s, and 

Lus{Nr±X"Y^) = {aj,Mr_r(gj){gju^xy)aj)>0, 
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where the latter inequaUty follows from Ms-r{g j){gj , xy) >: 0; see (|3.2|l . 
Similarly, one has 

i„=(iV;±r"Z'3) = iV>g ± ug„^>0, V|a + /3| <2r, 

because N^±Y°'Z'^ ~ Sifci '^k with degaK hx < 2s, and Ms_r{hK)i^K u^, yz) h 
0. Therefore, 

(3.8) Wi^ol < Nr y\a + f3\<2r; \u'oap\ < K V|« + /3| < 2r. 
In particular. 



inf Q, > 



NrY.\{Uy)o.f,0\+KY.\(fv^)0c.l3\ 

a(3 af3 



> —oo, 



provided s is sufficiently large. 

So, for sufficiently large s, let u'* be a nearly optimal solution of Q^, i.e., 

(3.9) inf < Lusif) < inf Q, + i < /* + i, 

s s 

and complete each sequence with zeros to make it an infinite sequence indexed in 
the canonical basis of M.[X, Y, Z]. Notice that by doing so, only elements of the form 
^^'-^ ''^Oa(i ^'"^ ^'-''^ zero. As (|3.8|l is true for arbitrary r, by a standard diagonal 
argument, there exists a subsequence {sfc} and an infinite sequence u indexed in 
the canonical basis of M.[X, Y, Z], such that 

(3-10) ^ "a/37 Va e N", (3 e N", 7 e N". 

From what precedes the only non zero elements Ua/37 of u, are those with a = or 
7 = 0. Next, introduce the subsequences 

(3.11) u'^y := {uaf}^ : 7 = 0}; := {u„;37 : a = }. 

Recall that the matrix Mr{u,xy) (resp. Mr{u,yz)) is identical to the moment 
matrix Mr{u^'^) (resp. Mr(u^^)) of the sequence u^^ (resp. w**^) indexed in the 
canonical basis of E[X, F] (resp. M[y, Z]). 

Similarly, as g,j € M.[X,Y] and /i^ G M[y, Z], the matrix Mr{gj u,xy) (resp. 
Mr{hK u,yz)) is identical to the localizing matrix Mrigju^^) (resp. Mr{hK u^^)) 
of the sequence m^*' (resp. u^^) indexed in the canonical basis of ]R[X, F] (resp. 
R[Y,Z]). 

Next, let r be fixed arbitrary. Then from the feasibility of u'^'' in Qs^, and 
the convergence (|3.1U|) . we obtain Mr{u,xy) = Mr{u^'^) >z and Mr(u,yz) = 
Mriu^^) >: 0. With same arguments, we also have 

Mr{gju,xy) = Af,(gj u"^) ^ 0, VJ C I,^ 
Mr{hKU,yz) ^ MrikKU^'') hO, yKCly^. 

As r was arbitrary, by Schmiidgen's Positivestellensatz fi}., it follows that the 
sequence u^*' (resp. u^^) has a representing measure /z^^ (resp. /i**^) with support 
contained in the compact set Vixy (resp. Ky^). Observe that 

u^o'p = "0/30 = ull, V/3 e N". 

Therefore, as measures on compact sets are moment determinate, /i^*' and /i**^ have 
same marginal on R™. 
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Next, the probability measure fi^^ on the cartesian product of Borel spaces 
M" X R™ can be disintegrated into a stochastic kernel q{dX \ Y) on M" given M™, 
and its marginal fi^ on M™, i.e., 



^,^yiAxB) = [ q{A\Y)iiy{dY), 

J B 



for all Borel rectangles [A x B) of R" x R™; see e.g. Bertsekas and Schreve |21 p. 
139-141]. Similarly, the probability measure /i^^ on the cartesian product of Borel 
spaces R™ x R^ can be disintegrated into a stochastic kernel q'{dZ \ Y) on R^* given 
R™, and its marginal on R™, i.e., 

^,y^iBxA) = I q'{A\Y)^jiy{dYi 



for all Borel rectangles {B x A) of R™ x R^. 

Let /J, be the probability measure on JR^+^+p defined by: 

(3.12) ^(AxBxC)= / q{A\Y)q\C\Y)^iy{dY), 

Jb 

for aU Borel rectangles (A x B x C) of M" x R™ x W. 
Taking C := Rp in (|TTa yields that 

fi{AxBxW)= I q{A\Y)^iy{dY) = ^r'y{Ax B), 

JB 

i.e., ^i^y is the marginal of /i on M"+™ (and in fact on Kj.j, because /.t^^(K^j^) = 1). 

Similarly, taking now A = M" , /x^^ is the marginal of /i on R™+p (and in fact on 
Kj^^). This clearly implies that /i is supported on K, i.e., fJ.(K.) = 1. Indeed, 

K = { (x, y, z) : {x, y) £ K^^y ; (y, z) e K^^ }, 

and so, 

m(K) - y q({x : {x,y) G K,,} | F) (7'({z : (y, z) £ K, J | F) /x^(dy) = 1, 

because from the definitions of /i^, K^-j, and Kj^^, 

<z({x : {x,y) E K,y}\Y) = q'{{z : (y, z) e K,,} | F) = 1 /x^-a.e. 
Finally, observe that from the convergence (|3.1U|) . we obtain 

(3.13) lim i„..(/) = Luif) = //d/.*, 
because / e R[X, Y] + R[Y, Z]. Therefore, by jn^J) and (E^SJ, we get 

/* > lim infQ^^, = lim L^'kif) = LuU) = f d^i. 

k — >oo k — ^oo J 

In view of H3.6|) and as ^ is supported on K, it follows that f* — J f dfi. Therefore, 
inf Qsj. — > /*, and as the sequence {Qr} is monotone nondecreasing, we obtain 

infQ. rr- 

Step (a). To prove absence of a duality gap between and its dual Q*, let v 
be the uniform probability measure on K, and let u be its sequence of moments. 
As K has nonempty interior, it follows that Mr{gju) >- 0, for every J C I^-j,, and 
all r, and similarly, Mrikx u) >- 0, for every K C Ij,^, and all r. 
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But this implies that AIr{gj u,xy) >- for all J C l^y and all r, because 
Mr{gj u,xy) is a submatrix of Mr{gju). Similarly, Mr{hK u,yz) >- for all 
K C and all r, because Mr{hK u, yz) is a submatrix of Mr{hK u). 

This means that m is a strictly feasible solution for and so, Slater's condition 
is satisfied for (see e.g. This in turn implies that there is no duality gap 

between and its dual Q*, i.e. supQ* — inf for all r > 0, and Q* is even 
solvable (i.e. supQ* = maxQ*) if has finite value; for more details on duality 
for semidefinitc programs, see e.g. Vandenberghe and Boyd 

Step (Hi). So let / be strictly positive on K, and let /* > be its global 
minimum on K, i.e., / > /* on K. From (i)-(ii), there exists some r such that 
maxQ* = inf Q,, > /*/2 > 0. Therefore, let A be an optimal solution of Q*. We 
have / - A e P{g) + P{h). But then / e P{g) + P{h) because A > 0, the desired 
result. 

(b) The proof of (b) is the same as that of (a), except that we now invoke 
Putinar's Positivstellensatz rather than Schmiidgen's. Indeed, if TV — ||(X, y)||^ £ 
Q{g): then Putinar's Positivstellensatz J holds, i.e., every polynomial of R[X, F], 
(strictly) positive on K^^j,, belongs to Q{g)\ see also Jacobi and Prestel ^.3';,. And 
so, the polynomial Nr ± X°'Y^ ^ strictly positive on J^xy for sufficiently large Nr^ 
belongs to the quadratic module Q{g) (instead of the preordering P{g)). The rest 
of the proof is identical. 

(c) The proof of (c) resembles that of (a), but the optimization problem is 
now the linear programming problem 

min £„(/) 

s.t. i„(.9"(l-g)/') >0, V(a,/3)eG, 

Lu{h°'il - hf) > 0, V (a, (3) eHr ' 
Mo = 1- 
where: 

Gr = {{a,P) e nI'-"I X NI'==«I : ^ deg^"^ +deggf < 2r} 

jeixy 

Hr = {(a,/3) e N''"-! X nI^"-! : ^ deg /i;!*- + deg /if'-' < 2r}. 

feeiyz 

The dual of the linear program is the linear program 

(3.15) L*: max{A : /-A e a(.g,l-g) + Cr(/i,l-/i)}, 

where Cr{g, 1 — g) C C{g, 1 — g) (resp. Cr{h, 1 — h) C C{h, 1 — hj) is the subcone 
Cr{g,l-g) = {aem[X,Y] : a ^ ^ c^^g'^il-gf; c„;5 > } 

{a,f3)GGr 

Cr{h,l-h)) = {a- e M[y,Z] : cr = ^ Ccph°'{l-hf- Ca/3>0}. 

By compactness of 'K.xy and Kj^^, the polynomial Nr^X°'Y^ is strictly positive on 
'Kxy, for some Nr^ and all a, (3 with |a + /3| < 2r. Therefore, using Krivine and 
Vasilescu JHl instead of Schmiidgen Positivstellensatz, there is some l{r) such that 
Nr±X°'YP e Ci{r){gA-9)- Proceeding as in (a) , the latter property is used in an 
optimal solution u'^ of L^, to bound Im^^qI ^^"^ + /^l < 2r, uniformly in s. We 

then obtain the convergence H3.1()|l for a subsequence {sk}- Using again Krivine 



(3.14) 
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and Vasilescu [ni], the subsequence of the hmit sequence u (see H3.11|l '). is the 
moment vector of a measure //^^ supported on 'i^xy^ and u^^ is the moment vector 
of a measure /i*'^ supported on Kj,^. In the present case, from linear programming 
duality, there is no duality gap between the linear programs and L*. The rest 
of the proof is along the same lines. 



3.3. Proof of Theorem I2.2L As Ka; and are both compact, then (possibly 
after some change of variables), we may and will assume that K^; C [—1/2, 1/2]" 
and Kj,z C [—1/2, 1/2]'"+^'. We only prove (a) because similar arguments hold for 
(b). Introduce the optimization problem: 

Luif) 



(3.16) 



mm 

u 

s.t. 



Mr{u,xy) >z 

^■^-^(gj) id J u,x) hO, V J C 



Mr 
Uq - 



r(hk){hK u, yz)>Q, y K C ly 



1. 



where now Mr{gj u, x) is the obvious analogue of Mr{gj u, xy) defined in tl3.1l The 
dual Q* now reads 



(3.17) q; 



max 

{o-j},{Vik},K0:A 

S.t. 



A 

/ 



K0 



ifCI„ 



e {aj} c S^, {^Pk} c 



degK0, deg o-jgj, degV'K^K < 2r 



Equivalently, Q* reads 

(3.18) Q; : max { A : / - A G (S^^), + Pr{g) + Pr{h) }, 

where (Sxy)r is the set of elements of S^y of degree at most 2r. 

Proceed as in the proof of Theorem l2.ir a). and so, let — (w^^^) be an arbitrary 
feasible solution of in (jXTB|) . As Ky, C [-1/2, 1/2]"+^, 1 ± V^Zl^ e P{h) for 
ah {a,P) e X W. Similarly, As K^^ C [-1/2,1/2]", 1 ± X" e P{g) for aU 
a £ N". Therefore, let r be fixed, arbitrary. With same arguments as in the proof 
of Theorem l2.1f a'). for all sufRciently large s. 



(3.19) 
(3.20) 



l<ool<l VaeN"; |a| < 2r 



< 1 



V(/3,7) e N" X NP, |/3 + 7|<2r. 



There is some additional technicality because we also need the boundedness of 



, uniformly in s. But from Ms{u,xy) > 0, we get u 



(2a)00 ^0(2/3)0 ^ 



«/3o)^ 



and so, in view of (|3.19|) - (|3.2U|) . 1 > |u^^q| for all s sufficiently large. Therefore, as 



in the proof of Theorem l2.1f a) let {u'*} be a sequence of nearly optimal solutions of 
Qg. The limit sequences u^^ and u^^ in (|3.11(l satisfy [w^al < 1 and 



\u%\ < 1 for all 



a, /3, 7. But this implies that u^^ is the moment sequence of a probability measure 
/Lt^^ supported on [—1,1]"+™; see Berg P Theor. 9]. In addition, again as in 
the proof of Theorem 12. If a) , and u^^ are moment sequences of two probability 
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measures fj,^ and with support contained in K.x and K.yz respectively. By 
construction, 

(3.21) j X'^ dii^ ^ ul^^ = u'^y^^ = j X'^ dii^y , VaeN". 

As both and /i^ have their support in a compact set, they are moment deter- 
minate, and so (|3.21|l impUes that /i^ is the marginal of fi^^ on K^;. Similarly, the 
marginal /i** of fjA^ on M™ is the same as the marginal of /i^^ on M'". Therefore, 
the measure defined in H3.12|l has marginal /i^ on K^;, marginal /i^^ on Kj,^, and 
marginal fi^^ on M" x R™. The rest of the proof is the same as that of Theorem 

inia). 

(b) One proceeds exactly as in (a), except that now one invokes Putinar's instead 
of Schmiidgen's Positivstellensatz. 
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